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ABSTRACT 

APEX single-dish observations at sub-millimeter wavelengths toward a sample of massive star-forming 
regions reveal that C2H is almost omni-present toward all covered evolutionary stages from Infrared Dark 
Clouds via High-Mass Protostellar Objects to Ultracompact Hn regions. High-resolution data from the 
Submillimeter Array toward one hot-core like High-Mass Protostellar Object show a shell-like distribution 
of C2H with a radius of ^9000 AU around the central submm peak position. These observed features 
are well reproduced by a ID cloud model with power-law density and temperature distributions and a 
gas-grain chemical network. The reactive C2H radical (ethynyl) is abundant from the onset of massive 
star formation, but later it is rapidly transformed to other molecules in the core center. In the outer 
cloud regions the abundance of C2H remains high due to constant replenishment of elemental carbon 
from CO being dissociated by the interstellar UV photons. We suggest that C2H may be a molecule well 
suited to study the initial conditions of massive star formation. 

Subject headings: stars: formation - astrochemistry - stars: individual (IRAS18089-1732) 



1. INTRODUCTION 

Spectral line surveys have revealed that high-mass star- 
forming regions are rich reservoirs of molecules from sim- 
ple diatomic s pecies t o complex and l arger molecules (e.g., 
Schilke et all Il997t lHatchell et all 119981: IComito et all 



20051 iBisschop et all 120071) . However, there have been 



rarely studies undertaken to investigate the chemical 
evolution during massive star formation from the ear- 
liest evolutionary stages, i.e., from High-Mass Starless 
Cores (HMSCs) and High-Mass Cores with embedded 
low- to intermediate-mass protostars destined to become 
massive stars, via High-Mass Protostellar Objects (HM- 
POs) to the final stars that are a ble to produce Ul- 
tracompact Hn regions (UCHns, see lBeuther et alj|2007l 
for a recent description of the evolutionary sequence). 
The first two evolutionary stages are found within so- 
called Infrared Dark Clouds (IRDCs). While for low- 
mass stars the chemical evolution from early molecu- 
lar freeze-out to more evolved protostellar cores is well 
studied (e.g.. iBergin fc Langerl 119971: iDutrey et "all 119971: 



iPavlvuchenkov et all 120061: IjQrgensen et all l2007f T. it is 
far from clear whether similar evolutionary patterns are 
present during massive star formation. 

To better understand the chemical evolution of high- 
mass star-forming regions we initiated a program to inves- 
tigate the chemical properties from IRDCs to UCHns from 
an observational and theoretical perspective. We start 
with single-dish line surveys toward a large sample obtain- 
ing their basic characteristics, and then perform detailed 
studies of selected sources using interferometers on smaller 
scales. These observations are accompanied by theoretical 
modeling of the chemical processes. Long-term goals are 
the chemical characterization of the evolutionary sequence 
in massive star formation, the development of chemical 
clocks, and the identification of molecules as astrophysi- 
cal tools to study the physical processes during different 
evolutionary stages. Here, we present an initial study of 
the reactive radical ethynyl (C2H) combining single-dish 



and interferometer observations with chemical modeling. 
Although C2H was previously observed in low-mass cores 
and Photon Dominated Regions (e.g., iMillar &: Freeman! 
1 1984t Hansen et allll995h . so far it was not systematically 
investigated in the framework of high-mass star formation. 

2. OBSERVATIONS 

The 21 massive star- forming regions were observed 
with the Atacama Pathfinder Experiment (APEX) in the 
875 fim window in fall 2006. We observed 1 GHz from 
338 to 339 GHz and 1 GHz in the image sideband from 
349 to 350GHz. The spectral resolution was O.lkms -1 , 
but we smoothed the data to ~0.9kms _1 . The average 
system temperatures were around 200 K, each source had 
on-source integration times between 5 and 16 min. The 
data were converted to main-beam temperatures with for- 
ward and beam effic iencies of 0.97 and 0.73, respectively 
(|Belloche et alJl2006t ). The average la rms was 0.4 K. The 
main spectral features of interest are the C2H lines around 
349.4 GHz with upper level excitation energies E u /k of 
42 K (line blends of C 2 H(4 5 5 - 3 4 4) & C 2 H(4 5 4 - 3 4 3) at 
349.338GHz, and C 2 H(4 4j4 '- 3 3 , 3 ) & C 2 H(4 4 , 3 - 3 3 , 2 ) at 
349.399GHz). The beam size was - 18". 

The original Submillimeter Array (SMA) C 2 H data 
toward the HMPO 1 8089-1732 were first presented in 
iBeuther et all |2005b). There we used the compact and 
extended configurations resulting in good images for all 
spectral lines except of C 2 H. For this project, we re- worked 
on these data only using the compact configuration. Be- 
cause the C 2 H emission is distributed on larger scales (see 
fj3|), we were now able to derive a C 2 H image. The inte- 
gration range was from 32 to 35kms , and the achieved 
la rms of the C2H image was 45 mJy beam -1 . For m ore 
details on these observations see IBeuther et all (<2005bD . 

3. RESULTS 

The sources were selected to cover all evolutionary 
stages from IRDCs via HMPOs to UCHns. We derived 
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our target lis t from t he samples of Klein et all (12005 ): 
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Fon tani et all (|2005h : Iffill et all (|2005l) : 
(200<1)! Table Q] lists the observed sources, their coordi- 
nates, distances, luminosities and a first order classifica- 
tion into the evolutionary sub-groups IRDCs, HMPOs and 
UCHns based on the previously available data. Although 
this classification is only based on a limited set of data, 
here we are just interested in general evolutionary trends. 
Hence, the division into the three main classes is sufficient. 

Figure [1] presents sample spectra toward one source of 
each evolutionary group. While we see several CH3OH 
lines as well as SO2 and H2CS toward some of the HM- 
POs and UCHns but not toward the IRDCs, the surpris- 
ing result of this comparison is the presence of the C2H 
lines around 349.4 GHz toward all source types from young 
IRDCs via the HMPOs to evolved UCHns. Table Q] lists 
the peak brightness temperatures, the integrated intensi- 
ties and the FWHM line-widths of the C2H line blend at 
349.399 GHz. The separation of the two lines of 1.375 MHz 
already corresponds to a line-width of 1.2 kms -1 . We have 
three C 2 H non-detections (2 IRDCs and 1 HMPO), how- 
ever, with no clear trend with respect to the distances or 
the luminosities (the latter comparison is only possible for 
the HMPOs). While IRDCs are on average colder than 
more evolved sources, and have lower brightness temper- 
atures, the non-detections are more probable due to the 
relatively low sensitivity of the short observations (H2J). 
Hence, the data indicate that the C2H lines are detected 
independent of the evolutionary stage of the sources in con- 
trast to the situation with other molecules. When compar- 
ing the line- widths between the different sub-groups, one 
finds only a marginal difference between the IRDCs and 
the HMPOs (the average Av of the two groups are 2.8 
and 3.1 km s -1 ). However, the UCHns exhibit significantly 
broader line- widths with an average value of 5.5 kms . 

Intrigued by this finding, we wanted to understand the 
C2H spatial structure during the different evolutionary 
stages. Therefore, we went back to a dataset obtained with 
the Submillimeter Array toward the hypercompact Hil re- 
gion IRAS 1 8089-1732 with a muc h higher spatial resolu- 
tion of ~ 1" (|Beuther et al.ll2005bD . Albeit this hypercom- 
pact Hil region belongs to the class of HMPOs, it is already 
in a relatively evolved stage and has formed a hot c ore with 
a rich molecular spectrum. iBeuther et al.l (|2005bT ) showed 
the spectral detection of the C2H lines toward this source, 
but they did not present any spatially resolved images. 
To recover large-scale structure, we restricted the data to 
those from the compact SMA configuration (|J2J). With this 
refinement, we were able to produce a spatially resolved 
C2H map of the line blend at 349.338 GHz with an angular 
resolution of 2.9" x 1.4" (corresponding to an average lin- 
ear resolution of 7700 AU at the given distance of 3.6 kpc). 
Figured] presents the integrated C2H emission with a con- 
tour overlay of the 860 /im continuum source outlining the 
position of the massive protostar. In contrast to almost all 
other molecular lines that p eak along with the dust contin- 
uum (IBeuther et al.ll2005bD . the C 2 H emission surrounds 
the continuum peak in a shell-like fashion. 

4. DISCUSSION AND CONCLUSIONS 

To understand the observations, we conducted a sim- 
ple chemical modeling of massive star-forming regions. A 



ID cloud model with a mass of 1200 M©, an outer ra- 
dius of 0.36 pc and a power-law density profile (p oc r p 
with p = —1.5) is the initially assumed configuration. 
Three cases are studied: (1) a cold isothermal cloud with 
T = 10K, (2) T = 50 K, and (3) a warm model with a 
temperature profile T oc r q with q = —0.4 and a tem- 
perature at the outer radius of 44 K. The cloud is illumi- 
nated by the interstellar UV radiation field (IRSF, iDraind 
I1978D and by cosmic ray particles (CRP). The ISRF at- 
tenuation by single-sized 0.1 /an silicate grains at a given 
radi us is calculat e d in a plane-parallel geometry follow- 
ing Ivan Dishoeckl (|1988| ) . The CRP ionization rate is as - 
sumed to be 1.3 x 10" 17 s" 1 (ISpitzer fc Tomaskd Il968h . 



The gas-grain chemical model by IVasvunin et all (2008) 
with the desorption energies and surface reactions from 
iGarrod et al.l (|2007t ) is used. Gas-phase rea ction rates are 
taken from RATE 06 (|Woodall et al.ll2007D . ini tial abun- 
dance s, were adopted from the "low metal" set of lLee et al.l 
(119981) . 

Figure [3] presents the C2H abundances for the three 
models at two different time steps: (a) 100 yr, and (b) 
in a more evolved stage after 5 x 10 4 yr. The C2H abun- 
dance is high toward the core center right from the be- 
ginning of the evolut i on, similar to previous models (e.g., 
Mill ar fc Neiadlll985t iHerbst fc Leunelll986t iTurner et all 
Il999f) . During the evolution, the C 2 H abundance stays 
approximately constant at the outer core edges, whereas 
it decreases by more than three orders of magnitude in 
the center, except for the cold T = 10 K model. The 
C 2 H abundance profiles for all three models show similar 
behavior. 

The chemical evolution of ethynyl is determined by rela- 
tive removal rates of carbon and oxygen atoms or ions into 
molecules like CO, OH, H 2 0. Light ionized hydrocarbons 
CH+ (n=2..5) are quickly formed by radiative association 
of C + with H 2 and hydrogen addition reactions: C + — > 

The protonated methane reacts 



CH+ 



CH 



CH+. 



with electrons, CO, C, OH, and more complex species at 
later stage and forms methane. The CH4 molecules un- 
dergo reactive collisions with C + , producing C2H2 anc ^ 
C 2 H^. An alternative way to produce C 2 H^ is the dis- 
sociative recombination of CH^" into CH3 followed by re- 
actions with C + . Finally, C2H2 an d C2H3" dissociatively 
recombine into CH, C2H, and C2H 2 . The major removal 
for C 2 H is either the direct neutral-neutral reaction with O 
that forms CO, or the same reaction but with heavier car- 
bon chain ions that are formed from C2H by subsequent 
insertion of carbon. At later times, depletion and gas- 
phase reactions with more complex species may enter into 
this cycle. At the cloud edge the interstellar UV radiation 
instantaneously dissociates CO despite its self-shielding, 
re-enriching the gas with elemental carbon. 

The transformation of C 2 H into CO and other species 
proceeds efficiently in dense regions, in particular in the 
"warm" model where endothermic reactions result in rich 
molecular complexity of the gas (see Fig. In contrast, 
in the "cold" 10 K model gas-grain interactions and surface 
reactions become important. As a result, a large fraction 
of oxygen is locked in water ice that is hard to desorb 
(£?des ~ 5500 K), while half of the elemental carbon goes 
to volatile methane ice (Ed cs ~ 1300 K). Upon CRP heat- 
ing of dust grains, this leads to much higher gas-phase 
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abundance of C2H in the cloud core for the cold model 
compared to the warm model. The effect is not that strong 
for less dense regions at larger radii from the center. 

Since the C2H emission is anti-correlated with the 
dust continuum emission in the case of IRAS 18089-1732 
(Fig-HJ), we do not have the H 2 column densities to quan- 
titatively compare the abundance profiles of IRAS 18089- 
1732 with our model. However, data and model al- 
low a qualitative comparison of the spatial structures. 
Estimating an exact evolutionary time for IRAS 18089- 
1732 is hardly possible, but based on the strong molec- 
ular line emission, its high centra l gas temperatures an d 
the ob served outflow-disk system (jBeuther et al.|[2004bl lal. 
2005b), an approximate age of 5 x 10 4 yr appears reason- 
able. Although dynamical and chemical times are not 
necessarily exactly the same, in high-mass star formation 
they should not diffe r t o much: Following the models by 
iMcKee fc Tar] (|2003l ) or lKrumholz et alj (|2007t) . the lumi- 
nosity rises strongly right from the onset of collapse which 
can be considered as a starting point for the chemical evo- 
lution. At the same time disks and outflows evolve, which 
should hence have similar time-scales. The diameter of 
the shell-like C 2 H structure in IRAS 18089-1732 is - 5" 
(Fig. [2]), or ^9000 AU in radius at the given distance of 
3.6 kpc. This value is well matched by the modeled re- 
gion with decreased C2H abundance (Fig. [3]). Although 
in principle optical depths and/or excitation effects could 
mimic the C2H morphology, we consider this as unlikely 
because the other observed molecules with many different 
transitions all peak toward t he central submm co ntinuum 
emission in IRAS 18089-1732 dBeuther et al.ll2005bh . Since 
C2H is the only exception in that rich dataset, chemical 



effects appear the more plausible explanation. 

The fact that we see C2H at the earliest and the later 
evolutionary stages can be explained by the reactive nature 
of C2H: it is produced quickly early on and gets replenished 
at the core edges by the UV photodissociation of CO. The 
inner "chemical" hole observed toward IRAS 18089-1732 
can be explained by C2H being consumed in the chemi- 
cal network forming CO and more complex molecules like 
larger carbon- hydrogen complexes and/or depletion. 

The data show that C2H is not suited to investigate 
the central gas cores in more evolved sources, however, 
our analysis indicates that C2H may be a suitable tracer 



N 2 H+ or NH 3 (e.g., Berrin et alj|2002t 


Tafalla et alJl2004t 


iBeuther et al.ll2005al iPillai et all 12006 


). While a spatial 



analysis of the line emission will give insights into the 
kinematics of the gas and also the evolutionary stage from 
chemical models, multiple C2H lines will even allow a tem- 
perature characterization. With its lowest J = 1 — 
transitions around 87 GHz, C2H has easily accessible spec- 
tral lines in several bands between the 3 mm and 850 /xm. 
Furthermore, even the 349 GHz lines presented here have 
still relatively low upper level excitation energies {E u /k ~ 
42 K), hence allowing to study cold cores even at sub- 
millimeter wavelengths. This prediction can further be 
proved via high spectral and spatial resolution observa- 
tions of different C2H lines toward young IRDCs. 
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Table 1 
Source parameters 



Name 


Type 


R.A. 


Dec. 


d 


L 




J T mb dv 


Atj(C 2 H) 
kms -1 


Ref. 






J2000 


J2000 


kpc 


log(L Q ) 


K 


Kkms" 1 




IRAS07029 


IRDC 


07:05:11.1 


-12:19:02 


1.0 


b 


0.21 


0.52 


2.4 ±0.3 


(1) 


IRAS08477 


IRDC 


08:49:32.9 


-44:10:47 


1.8 


b 


0.23 


0.76 


3.1 ± 0.5 


(2 3) 


TR AS09014 


IRDC 


09-03-09 8 


-47:48:28 


1.3 


b 








(2 3") 


TP A Q1 QOQQ 




1 0. 07-07 O 


£1 -9/1 -/1 7 
-01.24.4 ( 


9 1 

2.4 


b 










IK As 14000 


1KJJO 


14:03:36.0 


c 1 1000 

-61:18:28 


5.6 


b 


0.38 


1.13 


2.8 ± 0.5 


(2,3) 


TT) A OA0O1 1 


HlVtrO 


c\o 00 cro 


-42:07:57 


1.7 


3.5 


0.38 


0.85 


2.1 ± 0.2 


(2,3) 


TT) A OflO /I7H 

1KAS08470 


HlVtrO 


i~\0 AO AT CI 

08:48:47.9 


-42:54:22 


2.2 


4.2 


1.85 


6.81 


3.5 ± 0.4 


(3) 


TP A OnOKCO 

lKAb08563 


ilMrU 


no . c 0. 1 k 
08:58:12.0 


/i 0. 07.0 /i 
-42:3 /:34 


1 7 
1. f 


3.2 


1 no 
1.08 


o on 
3.30 


on 1 n -I 
2.9 ± 0. 1 


(2,3) 


TP A QAOl Q1 


TT1\/TPM 


uy. 14. 00.0 


-4f\o0.1o 


1 7 
1. 1 


o.4 


n 9Q 


O.OU 


i 10°. 
2.1 ± U . 




TR A ^0.090.0 

nx-riouyzuy 






-Ol.OO.ZO 


0.4 


A 1 










TP A QPIQc;70. 


TTA/TPO 


uy.oy.oi.u 


-0 1 .UO.40 


1 7 


. y 


n 91 
u.21 


U.04 


9 8 -t- fl A 
2.0 It U.4 




IRAS10123 


HMPO 


10:14:08.8 


-57:42:12 


0.9/3.0 c 


3.4/4.4 c 


0.17 a 


0.46 


2.5 ± 0.4 


(2 3) 


IRAS10184 


HMPO 


10:20:14.7 


-58:03:38 


5.4 


5.5 


0.41 


1.68 


3.9 ±0.3 


(3) 


IRAS 10276 


HMPO 


10:29:30.1 


-57:26:40 


5.9 


4.9 


0.24 a 


0.67 


2.6 ±0.5 


(3) 


IRAS 10295 


HMPO 


10:31:28.3 


-58:02:07 


5.0 


5.8 


0.69 


3.45 


4.7 ±0.3 


(3) 


IRAS 10320 


HMPO 


10:33:56.4 


-59:43:53 


9.1 


5.4 


0.85 


3.48 


3.8 ±0.6 


(3) 


G294.97 


UCH11 


11:39:09.0 


-63:28:38 


1.3/5.8 C 


3.9/5.3 c 


0.31 


0.78 


2.3 ±0.4 


(4) 


G305.20 


UCH11 


13:11:12.3 


-62:44:57 


3.0/6.8 c 


5.1/6.1 c 


0.44 


5.21 


11.1 ±0.5 


(4) 


G305.37 


UCH11 


13:12:36.3 


-62:33:39 


3.0/6.8 c 


d 


1.24 


6.28 


4.8 ±0.2 


(4) 


G305.561 


UCHn 


13:14:25.8 


-62:44:32 


4.0 


5.1 


0.94 


4.71 


4.7 ±0.6 


(4,5) 


IRAS14416 


UCH11 


14:45:22.0 


-59:49:39 


2.8 


5.1 


1.24 


6.19 


4.7 ±0.6 


(6) 



Ref.: (1) [Klein et all IpOol ). (2) iFontani et all (|2005l ). (3) lBeltran eTaTl (|2006f ). (4l lHill et all (|2005l ). (5) iFaundez et al l (|2004l ). (6) 
IVig et all (|2007l 1. 

a For these sources we list the parameters of the C2H line blend at 349.338 GHz, for all other sources it is the C2H line at 349.399 GHz. 
b Since the IRDCs are per default not detected at short wavelengths, they are no IRAS source and we cannot derive a luminosity. 
c Near and far distances and corresponding luminosities. 
d No IRAS counterpart, hence no luminosity estimate. 
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Fig. 1. — Sample spectra obtained with APEX in a double-sideband mode. The spectra cover 1 GHz of data around 349.5 GHz and 1 GHz 
around 338.5 GHz. All spectral lines are labeled, lines in the lower sideband are marked with a "(L)". The top-panel shows an IRDC example, 
the middle-panel a typical HMPO/hot core and the bottom panel a UCHll region. The C2H line-widths are indicated in each panel. 



Fig. 2. — The grey-scale shows the integrated emission (from 32 to SSkms" 1 ) of the line blend of C 2 H(4 5 , 5 - 34.4) and C 2 H(4 5 , 4 - 3 4 ,3) 
around 349.338 GHz obtained with the Submillimeter Array using only the compact configuration data (Bcuthcr ct al. 2005b). The resulting 
synthesized beam is shown at the bottom left corner (2.9" X 1.4"). The C2H emission is presented in grey-scale with thin contours (dashed 
contours negative features), and the 860 fim continuum peak is shown in thick contours. The C2H emission starts at the 2<r level and 
continues in lcr steps with the la level of 450 m Jy beam - . The 860 /im emission is contoured from 10 to 90% (step 10%) of the peak emission 
of 1.4Jybeam _1 . 
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Fig. 3. — Top left: Radial profiles of the C2H abundance for the three cloud models at early times, t < 100 yr. Bottom left: The same but 
for the later time step of 5 X 10 4 yr. The parameter q denotes the exponent of the assumed temperature distribution (see also main text). 
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